Further Evidence for Neutrino Flux Variability from Super-Kamiokande Data 
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While KamLAND apparently rules out Resonant-Spin-Flavor-Precession (RSFP) as an expla- 
nation of the solar neutrino deficit, the solar neutrino fluxes in the CI and Ga experiments vary 
with solar rotation rates. Added to this evidence, summarized here, a power spectrum analysis of 
the Super-Kamiokande data reveals (99.9% CL) an oscillation in the band of twice the equatorial 
rotation frequencies of the solar interior. An m — 2 magnetic structure and RSFP, perhaps as a 
subdominant process, would give this effect. Solar cycle data changes are seen, as expected for 
convection zone modulations. 

PACS numbers: 26.65.+t, 95.75.Wx, 14.60.St 



Recent results from the KamLAND experiment |l| 
seem to confirm the Large-Mixing Angle (LMA) solution 
to the solar neutrino deficit and rule out the Resonant- 
Spin-Flavor-Precession (RSFP) explanation 0. On the 
other hand, there is increasing evidence 0-0 that the so- 
lar neutrino flux is not constant as assumed for the LMA 
solution, but rather it varies with well known solar rota- 
tion periods. The solar neutrino situation may be com- 
plex, and Spin-Flavor-Precession could be subdominant 
to LMA, as suggested in 0, or if there is at least one 
light sterile neutrino even RSFP could be a subdominant 
process. Since this recent information on solar neutrino 
variability is not widely known, a brief summary is pre- 
sented of analyses of radiochemical neutrino data, along 
with new input from the Super-Kamiokande experiment 
. While the 10-day averages of Super-Kamiokande so- 
lar neutrino data ^| show no obvious time dependence, 
a power-spectrum analysis |l 3| displays a strong peak at 
the frequency 26.57 ± 0.05 y _1 (period 13.75 d), where 
the width of peaks are computed for a probability drop to 
10%. The probability of finding this peak (or a stronger 
peak) by chance within a band specified by twice the 
near-equatorial rotation frequencies of the solar interior 
(25.36-27.66 y" 1 ) is found to be 0.001 (99.9% CL). This 
frequency is also seen in the radiochemical data. 

Whether or not it is the correct interpretation, the 
RSFP framework provides a simple way to understand 
the data, and hence it will be used here. The RSFP mech- 
anism, which requires a neutrino transition magnetic mo- 
ment, provides at least as good a global fit jl4lll5l | to solar 
data (which depends mainly on the Super-Kamiokande 
spectrum) and a better fit to average rates of the indi- 
vidual experiments than does LMA. This is because the 
neutrino survival probability, while having a resonance 
pit at a density that suppresses the 0.86 MeV 7 Be line 
(as does the Small-Mixing- Angle (SMA) solution), tends 
at high energies toward 1/2 and hence fits the Super- 
Kamiokande spectrum, whereas the survival probability 



goes to unity in the SMA case. Thus, as a subdominant 
process, it could also improve fits to the data. 

A brief review follows of the published evidence for so- 
lar neutrino flux variability, put into a coherent neutrino 
scheme. Byanalyzing 10 3 simulated data sequences, it 
was found [3| that the variance of the Homestake solar 
neutrino data is larger than expected at the 99.9% 
CL. A power spectrum analysis of the data showed a 
peak at 12.88 ± 0.02 y" 1 (28.4 d), compatible with the 
rotation rate of the solar radiative zone. Four sidebands 
gave evidence at the 99.8% CL for a latitudinal effect as- 
sociated with the tilt of the sun's rotation axis, and the 
latitude dependence was also seen directly in the data at 
the 98% CL 0. The GALLEX data [3 showed a peak at 
13.59 ±0.04 y" 1 , compatible with the equatorial rotation 
rate of the deep convection zone. This peak is also in the 
Homestake data, and a combined analysis of both data 
sets shows that the 13.59 y" 1 peak is larger than in either 
data set alone. Comparison [7J of the power spectrum for 
the GALLEX data with a probability distribution func- 
tion for the synodic rotation frequency as a function of 
radius and latitude, derived from SOHO/MDI hclioseis- 
mology data [l^, results in Fig.^ This map shows the 
rotation frequency coincides with the neutrino modula- 
tion in the equatorial section of the convection zone at 
about 0.8 of the solar radius, Rq. The influence of these 
rotation frequencies extends even to the corona, since 
the SXT instrument on the Yohkoh spacecraft pro- 
vides X-ray evidence for two "rigid" rotation rates, one 
(13.55 ± 0.02 y _1 ) mainly at the equator, and the other 
(12.86±0.02 y _1 ) mainly at high latitudes. These values 
are in remarkable agreement [8j with the neutrino modu- 
lation frequencies and also with their equatorial location 
in one case and non-equatorial in the other; see Fig. [5] 

That two dominant frequencies are shown by both 
coronal X-rays and neutrino flux is a feature of magnetic 
rotation well known in solar physics. For instance, an 
analysis of the photospheric magnetic field during so- 
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FIG. 1: Map of the resonance statistic of the SOHO/MDI 
helioseismology and GALLEX data on a meridional section 
of the solar interior. The only high probability areas (red) 
are lens-shaped sections near the equator, and all others are 
low probability. 
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FIG. 2: Comparison of normalized probability distribution 
functions formed from power spectra of data from SXT equa- 
tor (red), SXT N60-S60 (green), Homestake (black), and 
GALLEX (blue). Note that the SXT (red) and GALLEX 
data are equatorial, and the other two are not. 



lar cycle 21 found two dominant magnetic regions: one in 
the northern hemisphere with synodic rotation frequency 
~ 13.6 y _1 , and the other in the southern hemisphere 
with synodic rotation frequency ~ 13.0 y _1 . Similarly, 
an analysis (2l| of flares during solar cycle 23 found a 
dominant synodic frequency of ~ 13.5 y _1 for the north- 
ern hemisphere and ~ 12.9 y _1 for the southern hemi- 
sphere. These and other studies show a strong tendency 
for magnetic structure to rotate either at about 12.9 y _1 
or at about 13.6 y . Since the sun's magnetic field is 



believed to originate in a dynamo process at or near the 
tachocline, it is possible that some magnetic flux is an- 
chored in the radiative zone just below the tachocline, 
where the synodic rotation frequency is about 12.9 y , 
and some just above the tachocline in the convection 
zone, where the synodic rotation frequency is about 13.6 
y -1 . It is also possible that the 12.9 y _1 frequency re- 
sults from a latitudinal wave motion in the convection 
zone excited by structures at or near the tachocline. 

An example of latitudinal oscillatory motion of mag- 
netic structures may be the well-known Rieger-type os- 
cillations with frequencies of about 2.4, 4.7, and 7.1 
y _1 [23. These periodicities may be attributed to r- 
mode oscillations with spherical harmonic indices £ = 3, 
m = 1,2,3, giving frequencies v = 2mvn/£{£ + 1) which 
are seen in neutrino data 0.0. Q A joint spectrum anal- 
ysis of Homestake and GALLEX-GNO data yields 
peaks at 12.88, 2.33, 4.62, and 6.94 y -1 , indicating a 
sidereal rotation frequency vr — 13.88 ±0.03 y _1 . While 
I > 2 is required, odd— I values have nonzero poloidal ve- 
locity at the equator that could move magnetic regions 
in or out of the neutrino beam to earth. 

The difference between the main modulations detected 
by the CI and Ga experiments may be explained by 
the tilt of the solar axis relative to the ecliptic, along 
with the fact that CI and Ga neutrinos are produced 
mainly at quite different radii. The Ga data comes 
mostly — especially as the fit requires suppressing the 7 Be 
line — from pp neutrinos, which originate predominantly 
at large solar radius (~ 0.2 R©), so that the wide beam 
of neutrinos detected on earth is insensitive to axis tilt. 
Thus the beam of neutrinos detected by the Ga experi- 
ments exhibits no seasonal variation and can be modu- 
lated by the equatorial structure indicated in Fig.^ lead- 
ing to the observed frequency at about 13.6 y -1 . On the 
other hand, the Homestake experiment detects neutrinos 
produced from a smaller sphere (~ 0.05 R©), so that the 
axis tilt causes these neutrinos mainly to miss the equa- 
torial structure of Fig.^and instead sample nonzero lat- 
itudes where the 12.9 y _1 modulation may occur. Twice 
a year axis tilt has no effect for these neutrinos, leading 
to a seasonal variation in the measured flux Q . 

The 13.6 y _1 frequency, located as in Fig.^ represents 
a modulation of the pp neutrinos which are at or near the 
steeply falling edge of the neutrino survival probability. 
This is close to the RSFP resonance pit suppressing the 
7 Be neutrinos, which is where the largest value of Am 2 /E 
satisfies 



Am 2 /E = 2V2G F N< 



off, 



(1) 



and is essentially the same as for an MSW resonance 
For MSW, ATeff = N e (the electron density), and 
for RSFP N c g = N e — N n for Majorana neutrinos, where 
N n is the neutron density (about N e /6 in the region of 
interest in the Sun). For Dirac neutrinos N e g = N e — 
N n /2, but only Majorana neutrinos provide a fit to solar 
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FIG. 3: Cumulative Rayleigh power vs. time for the 13.59 y 
frequency peak from the GALLEX-GNO data. Note that the 
power builds up from the start of data taking after the 1990 
solar maximum until the 1996.8 solar minimum, after which 
there is little evidence for that frequency. 
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FIG. 4: Normalized GALLEX flux measurement runs (prior 
to 1997) reordered according to the phase of the 13.59 y _1 
solar modulation. The division in phase is made so as to 
have equal numbers of events in the descending (blue) and 
ascending (red) parts of the cycle. 



experimental rates with RSFP [lj, [15| , a very important 
consequence of such a solar solution. Using the N e and 
N n values at the 0.8 R© location of Fig.^in Eq.^results 
in Am 2 /E ~ 10~ 14 eV, putting the RSFP resonance pit 
at the location needed to fit the solar data. N e g varies 
exponentially with radius and would be quite different 
for a somewhat changed neutrino modulation frequency. 

It is well known that the convection zone magnetic field 
changes with the solar cycle, so the neutrino modulation 
features described above are not permanent. It has been 
argued |24| that an RSFP effect would have to be in the 
radiative zone, where the field does not change with the 
solar cycle, under the assumption solar cycle variations 
are not observed. On the contrary, we show here that 
solar cycle changes play an important role. Variations in 
neutrino rates are difficult to observe, since changes in 
field magnitude would be undetectable if the transition 
remains adiabatic, and even if flux modulation results, 
average rates may vary only slightly. The feature that 
is most sensitive to field magnitude or radial variations 
is the intersection of the very steeply falling pp neutrino 
spectrum with the falling edge of the RSFP resonance pit. 
As can be seen in Fig. |3 the 13.6 y^ 1 frequency associ- 
ated with these neutrinos increased in strength from the 
start of data taking after the solar maximum of 1989.6 
to the solar minimum of 1996.8, after which the modula- 
tion becomes weak. Also it was during that cycle that the 
main buildup in the strength of the 12.9 y" 1 oscillation 
was detected by Homestake. The SXT X-ray data, with 
which these two frequencies had remarkable agreement 
8] , also came from that same solar cycle. 
Another feature attributable to a time variation of the 



intersection of the pp spectrum with the edge of the RSFP 
resonance pit is the appearance of a bimodal flux distri- 
bution in the Ga data. This distribution is clearly evident 
during the same solar maximum to solar minimum 
period when the flux is binned appropriately using indi- 
vidual runs, instead of averages over several data runs. 
This neutrino flux effect also diminishes after the solar 
minimum. Adding to the significance of this result is the 
plot of Fig. 01 which shows that when the end times of 
runs are reordered according to the phase of the 13.59 
y" 1 rotation, the flux values are low in one-half of the 
cycle and high in the other. 

It is unfortunate that the Homestake experiment, 
which detected mainly the 8 B neutrinos (the only 
neutrino component registered by Super-Kamiokande) 
stopped operating at about the time that Super- 
Kamiokande started. As a result, there is no way to 
predict from results of other experiments what neu- 
trino flux variation should have been detected by Super- 
Kamiokande during its operation from May 1996 (near 
solar minimum) until July 2001 (near solar maximum). 

Recently the Super-Kamiokande group released [l^ 
flux measurements in 184 bins of about 10 days each. 
While these measurements vary by ~ 2, averaged over 
all bins their fractional error is 0.14. Because of the 
regularity of the binning, the "window spectrum" (the 
power spectrum of the acquisition times) has a huge peak 
(power S > 120) at a frequency of v = 35.98 y _1 (pe- 
riod 10.15 d). (Note that the probability of obtaining a 
power of strength S or more by chance at a specified fre- 
quency is e~ s .) This regularity in binning leads to alias- 
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FIG. 5: Combined power spectrum of the Super-Kamiokande, 
SAGE, and GALLEX-GNO data. 

ing of the power spectrum of the flux measurements. In a 
spectrum formed by a likelihood procedure, the strongest 
peak in the range to 100 y _1 occurs at v = 26.57 y _1 
[Hill with 5 = 11.26. In the range 0-40 y _1 , the next 
strongest peak is at v — 9.41 y" 1 with S = 7.33. Since 
26.57 + 9.41 = 35.98, we infer that the weaker peak (9.41 
y _1 ) is an alias of the stronger (26.57 y _1 ). 

Having found this frequency in the high-statistics 
Super-Kamiokande data, we then examined the power 
spectrum from a combined analysis of the SAGE [26| and 
GALLEX-GNO data, the sampling time of the Home- 
stake data being too long for this high a frequency. 
The combined experiments showed a peak at 26.54±0.05 
y _1 with S = 5.75. When combined with the Super- 
Kamiokande data, we find a peak at 26.54±0.03 y _1 with 
S = 13.95, as shown in Fig. There is no significant 
peak near 9.41 y _1 , further showing this to be an alias 
frequency. 

Since 26.5 y~ x is within the band of twice the near- 
equatorial rotation frequencies in the solar interior, this 
result indicates that some 8 B neutrinos are experiencing 
an "m = 2" structure, two circumferential regions of the 
magnetic field which differ in strength from the rest. This 
modulation, unlike the Gapp threshold effect, is not very 
deep (~ 10%). Such oscillations at the harmonic of the 
rotation frequency are not uncommon in solar data. 

The analysis of Super-Kamiokande data, when com- 
bined with the results of analyses of CI and Ga data, 
yields evidence for variability of the solar neutrino flux. 
The KamLAND experiment supports the LMA solution, 
but it could be combined with RSFP. A definitive result 
on neutrino flux time dependence can be carried out using 
a combined analysis from all experiments, particularly 
using one-day bins from all water Cerenkov detectors. 
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